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o co 

.p-  „.| — - For  alracet  two  centuries  many  quantitative  experiments  have  been  made 

" *-  concerning  the  flow  of  water  down  an  inclined  open  channel.  HydraulicUna 
have  evolved  empirical  laws  in  confusing  multiplicity  to  describe  a uniform 
flow  as  a function  of  the  slope  and  the  resistance  caused  by  friction  with  the 
channel  walls  and  by  turbulence.  Each  of  the  laws  has  its  own  best  rnrge  of 
applicability,  but  most  &r*>  of  the  same  basic  type,  vb.  the  effective  resistance 
is  considered  in  the  momentum  equation  as  an  opposing  body  force  of  form 
— Xu’/y",  u being  the  velocity,  y the  “hydraulic  radius”  and  X a coefficient 
describing  the  channel  roughness,  which  in  this  notation  is  independent  of  u 
and  y.  The  exponents  m,  n > 0 have  specific  values  characteristic  of  the  differ- 
ent formulas.  For  a wide  rectangular  channel,  y will  be  also  the  depth  of  the 
water. 


A uniform  (turbulent)  flow  with  a free  surface,  being  subject  to  some  such 
resistance  law,  sometimes  reaches  a condition  of  instability,  then  changes  into 
a complicated  progressing  wave  pattern  called  “roll  v,avc3”  which  are  periodic 
in  distance,  with  the  profile  moving  downstream  faster  than  the  water  particles. 
Mathematical  analyses  have  been  made  of  certain  aspects  of  the  final  roll  wave 
motion  by  Thomas  [1]  and  one  of  the  present  authors  [21.  Results  giving  neces- 
sary criteria  have  been  obtained  by  considering  the  correct  elope  o:  the  roll  wave 
profile  or  by  satisfying  the  energy  inequality  at  the  moving  discontinuities 
(bores).  Our  present  purpose,  however,  is  rather  to  investigate  the  instability 
of  the  original  uniform  flow  in  its  dependence  upon  the  specific  resistance  function 
assumed,  following  the  method  of  Jeffreys  [3]  who  first  did  this  in  1925  for  the 
Chezy  formula  (n  — 2,  m — 1 in  our  present  notation).  The  criterion  obtained 
by  Jeffreys  for  in?itKbility  is  identical  with  the  condition  later  obtained  for  roll 
wave  formation  in  [1]  and  [2].  Each  of  the  papers  [1],  [2]  and  [3]  used  the  Chezy 
resistance  fuuctiou. 


In  1940  Keulegan  and  Patterson  [4]  derived  a stability  criterion  for  the 
Manning  formula,  based  upon  an  expression  for  wave  celerity  due  to  Boussinesq. 
Another  type  of  stability  argument  was  later  used  by  Vodemikov  [51,  using 
certain  approximations  of  Saint  Venant.  His  method  is  rather  complicated, 
but  the  results  are  applicable  for  a general  resistance  term.  More  recently 
Craya  [6]  has  refined  this  approach  by  considering  the  growth  of  an  infinitesimal 
Bhock  (“elementary  wave”)  with  results  applying  also  in  the  case  of  a general 
resistance.  From  the  qualitative  standpoint,  Cornish  [7],  with  appendix  by 
Harold  Jeffreys,  has  presented  some  interesting  remarks  on  the  formation  of  roll 
waves  with  respect  to  the  resistive  action  of  the  stream  bed. 
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Using  the  shock  energy  inequality  to  construct  actual  roll  wave  eolutiooa, 
it  was  observed  in  [2]  that  such  solutions  could  never  be  constructed  if  the  varia- 
tion of  resistance  effect  with  depth  was  ignored.  That  is,  if  the  resistance  de- 
pended only  on  the  square  of  the  velocity,  the  energy  shock  condition  for  the 
bores  of  the  waves  would  always  be  violated.  Such  a simplifying  assumption 
on  the  resistance  io  sometimes  made,  however,  particularly  when  variations  in 
depth  are  considered  to  be  relatively  small.  We  wish  now  to  determine  whether 
agreement  between  the  stability  approach  and  the  roll  wave  approach  extends 
also  to  this  case,  and  to  see  more  generally  how  instability  depends  on  the  form 
of  the  resistance  function.  That  is,  one  wishes  to  aee  what  particular  quantities 
in  the  resistance  mechanism  can  lead  to  instability  and  roll  waves.  This  present 
investigation  will  adapt  the  method  of  Jeffreys  to  2 general  resistance  function, 
and  results  will  be  compared  with  those  obtained  by  the  other  authors  mentioned 
above,  using  different  approaches  to  tV  problem. 

In  a two-dimensional  flow,  let  u be  the  velocity  component  parallel  to  the 
stream  bed,  inclined  at  an  angle  $ below  the  horizontal,  and  let  y be  the  depth  of 
water.  The  non-linear  shallow  water  equations  are  then 

u,  + uu.  + y cos  0y.  = y Bin  0 - (Xu"/y“) 

(1) 

yu,  + uy,  -f  y,  - 0. 

The  Chexy  formula  corresponds  to  n - 1,  n « 2;  the  Manning  formula  to 
m — 4/3, » «■  2;  while  in  the  various  Lea  and  Barnes  expressions,  m varies  from 
1.13  to  1.64  and  n from  1.7  to  2.15  (see  King  [8]). 

The  poesible  uniform  flows  are  given  by  u — U,  y ■»  Y where 

» -cw- 


Considering  small  deviations  given  by  u — U + u(x,t),  y — Y + y(x,<)  and 
line*n ring  yield  two  equations  for  u and  y, 


gm  sin  8 _ ..  nlT~  _ 

u,  + Liu,  + Q cos  9y,  - * — p — y - A — — u 


(3) 


Uy.  + Yu,  + y,  - 0. 

When  y is  eliminated,  we  get  the  second  order  linear  equation  for  v(zf). 


u„  + 2UU.,  + (U*  - gY<x*ti^  + t2^-u, 


(4) 


( . . , XnLT-'N  _ _ 

T \mg  sm  I + — : ~ ) «.  - 0. 


We  now  study  the  time  growth  or  decay  of  the  class  of  progressing  wave  aoh>- 


RESISTANCE  EFTECTS  ON  HYDRAULIC  INSTABILITY 


95 


tions  of  type 


with  wave  speed  a/0  And  wave  length  proportion  a!  to  i/0.  Let  y “ r + is;  then 
insertion  of  the  boiutioua  A exp  [yl  + i&i  } into  (4)  implies  the  following  relation 

between  y and  0 

\nir~'  . /xVt/*-*  ^ „ ...  . „Y" 

(5)  y tU0  ic  P qY  cob  8 - 0gm  sin  8)  . 

For  stability  we  require  both  roots  to  be  in  the  Jef ' ha'i  of  the  complex  y~ 
plane.  If  parameters  are  adjusted  bo  that  the  real  part  of  the  first  square  root  in 

(5)  is  equal  to  X»tT'*/ 2K",  then  the  corresponding  y value  will  be  on  the  imagi- 

nary axis,  and  the  other  value  in  the  left  haif-piane.  We  therefore  equate  this 
square  root  to  XnLT~ 1 fit  -f-  it,  separate  reals  and  imaginarica,  eliminate  6,  and 

use  relation  (2).  The  condition  for  the  right  root  to  be  on  the  s-oxis  is  then 

(6)  mV*  yt~  ibt  f . n-x*  ^ 9 

When  the  expression  on  the  left  i?  greater  or  smaller,  it  can  be  shown  that  the 
root  moves  to  the  left  or  right,  respectively.  Considering  the  expressions  in 

(6)  as  functions  at  8,  we  see  that  equation  (6)  must  have  exactly  one  8 solution 
between  0 and  w/2.  The  stability  criterion  tor  general  resistance  is  thus 

(7)  aT  K*~  sin*  • | n**X*  cos"  8 

implying  stability  for  the  two  top  signs,  and  instability  for  the  bottom  sign. 
The  condition  for  instability  can  be  put  into  an  equivalent  form,  using  (2), 

(8)  U>„  VVFcos#, 

which  is  identical  with  the  Vtdemikov  and  Craya  results,  obtained  by  other 
methods  and  assumptions. 

The  above  instability  criterion  for  the  Cbesy  formula  reduces  to  tan  # > 
4X,  the  result  obtained  in  (1),  (2],  and  (3).  For  the  Manning  formula  (vt  » 4/3, 
w — 2),  our  result  is 

(9)  U*  > | gY  cos  8 

which  agrees  with  the  result  of  Keulegan  and  Patterson  except  for  the  absence 
of  the  factor  cos  8 in  their  criterion. 

To  see  what  is  tl»  essential  non-linearity  in  (1)  causing  instability,  we  now 
consider  two  q>edal  limit  eases  of  (7): 
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Cate  1 : m — ► 0,  n > 0. 

Now  (7)  reduces  to  only  one  possibility 

0 < n*"X*  cos*  6, 

hence  we  conclude  that  any  flow  governed  by  a resistance  varying  with  any 
power  of  velocity,  but  independent  of  depth,  must  remain  stable. 

Cate  2:  n — ► 0,  m > 0. 

By  virtue  of  the  law  for  uniform  regime  Y — (\/g  sin  0),/*J  the  only  possi- 
bility for  (7)  is  the  equality 

g * Ylm  sin*  0 •**  X*. 

Both  y loots  move  to  the  x-axis,  still  a region  of  stability.  Although  this  case, 
describing  a type  of  hydraulic  Coulomb  friction,  may  not  correspond  to  actual 
flows,  the  analysis  nevertheless  indicates  that  both  effects  of  depth  and  velocity 
dependence  operating  simultaneously  are  needed  to  produce  instability.  For 
this  reason,  since  one-dimensional  compressible  gas  flow  is  governed  by  the 
same  equations,  with  resistance  depending  only  on  velocity,  one  would  not 
expect  any  analogous  type  of  instability  to  develop  in  such  flows. 
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